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Recent PAMELA data show that positron fraction has an excess above several GeV 
while anti-proton one is not. Moreover ATIC data indicates that electron/positron 
flux have a bump from 300 GeV to 800 GeV. Both annihilating dark matter (DM) 
with large boost factor and decaying DM with the life around 10^^ s can account 
for the PAMELA and ATIC observations if their main final products are charged 
leptons (e, /i and r). In this work, we calculated the neutrino flux arising from /x and 
T which originate from annihilating/decaying DM, and estimated the flnal muon rate 
in the neutrino telescopes, namely Antares and IceCube. Given the excellent angular 
resolution, Antares and IceCube are promising to discover the neutrino signals from 
Galactic center and/or large DM subhalo in annihilating DM scenario, but very 
challenging in decaying DM scenario. 

PACS numbers: 95.35.+d,13.15.+g,95.55.Vj, 98.62.Gq 
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I. INTRODUCTION 

The main components of our Universe are 'dark'. Many astrophysical observations have 
confirmed the existence of dark matter (DM) which contributes roughly 23% to the energy 
density of the Universe. If the DM are weakly interacting massive particles (WIMPs), 
they would annihilate into gamma rays, neutrinos and anti-matter particles, which could be 
detected experimentally. Generally speaking, the DM should be stable compared with the 
life of the Universe. However it is possible that the DM is decaying at an extremely low 
rate. Similarly, the DM decay products can be observable. 

Recently PAMELA released their first cosmic rays results on the positron and anti-proton 
ratios Usually the anti- matter particles are expected to be produced when the cosmic 
rays propagate in the Galaxy and interact with the interstellar medium. The PAMELA 
data show a clear upturn (excess) on the positron ratio above ~ 10 GeV up to ~ 100 GeV. 
Besides PAMELA observation, ATIC collaboration also reported an electron/positron excess 
around 300 ~ 800 GeV |2|. Contrary to the excess of position ratio, PAMEMA anti-proton 
ratio is consistent with the predictions. Both PAMELA and ATIC observations indicate 
some unknown high energy primary positron sources in the Galaxy. Studies showed that 
the pulsars or other nearby astrophysical sources may account for the PAMELA results [sl . 

Other possibility to induce the positron excess can be the annihilating 4,[5, 6] or decaying 
{3] DM in the Galaxy. In the annihilating DM scenario, if the cross section (erf) is taken 
as the typical annihilation cross section ~ 3 x 10~^^ cm^s~^ for thermal relic, a large extra 
boost factor (BP) is needed to produce sufficient positron flux to fit the PAMELA results. 
Many ideas are proposed to induce the large BP. For example, the new 'long' range force 



would enhance the annihi 
'Sommerfeld effect 



ation cross section for today low velocity DM, which is dubbed as 



For the decaying DM scenario, by adjusting the lifetime of 
the DM, no BP is needed. Contrary to positron excess, no obvious anti-proton excess has 
been observed. This suggests that the main final products from the annihilation or decay 
of DM are charge leptons (e, fi and r) for m^A/ < ITeV. Otherwise one needs much heavy 
DM, say O{10TeV) 5|], or some special cosmic-ray propagation parameters. 

Provided that charged leptons are generated by DM annihilating and/or decaying DM, 
it is quite interesting to investigate which charged leptons are response for the PAMELA 
and ATIC observations. It is natural to expect that the high energy charged leptons are 
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correlated to other observable signals, such as gamma rays, synchrotron emissions or neu- 
trinos. Detecting neutrino signals from DM is very attractive and useful to pin down the 
nature of DM. In this paper, we will calculate the neutrino flux from the fi and r which are 
producing by DM. Unlike other DM related products, neutrinos have less trajectory defec- 
tion and less but finite energy loss 10| during propagation, due to their weakly interaction 
nature. As a result, the neutrinos can keep the important information of the DM. For the 
same reason, the neutrinos are difficult to be detected if the flux is low (for some recent 
works, see Inl, Q, Q, 14 , 15 



such as Super-Kamiokande 



16| ). Usually the neutrinos detected by neutrino telescopes, 
171 . Il8| , AMANDA 19| etc. are expected to be the atmospheric 



neutrinos, which are produced by interaction between the cosmic rays and atmosphere. The 
atmospheric neutrinos are the main backgrounds for detecting the extra neutrinos from DM. 
Currently, no high energy astrophysical neutrinos have ever been detected 2o| . 

If the charged leptons are indeed produced by DM, it is a good news for neutrino detection. 
Firstly, the neutrino flux can be much larger than that of usual expected. The PAMELA 
and ATIC data need higher primary positron production rate than usual production, thus it 
is natural to expect that large neutrino flux can be produced from /i and r decays. Secondly, 
the heavy DM (e.g. ~ TeV) is favored by ATIC result, thus the neutrinos at the similar 
energy scale are expected. As the main backgrounds, the atmospheric neutrino flux decreases 
logarithmically with the increment of energy, thus higher energy neutrinos correspond to less 
backgrounds. On-going and future kilometer size neutrino telescopes such as Antares 21| . 

can better explore these signals. Thirdly, in the annihilating 



IceCube 



22[ and KMSNeT 



DM scenario, the flux of neutrinos depends on the square of the DM number density. If 
one searches the Galactic center (GC) or large DM subhalos which contain the denser DM 
than other areas, large neutrino flux can be produced. Moreover if the Sommerfeld effect 
of DM annihilation is true, the DM annihilation in the subhalo can acquire even larger BF 
due to the lower velocity dispersion than that in the DM halo . As a side remark, this 
feature provides a potential method to distinguish two DM scenarios, namely GC or DM 
subhalos can be higher flux neutrino sources in the annihilating DM scenario than those in 
the decaying DM one. 

This paper is organized as following. In section II, we calculated the neutrino flux from 
the GC and DM subhalo, in the model-independent manner, in the annihilating and decaying 
DM scenarios. In section III, we estimated the muon rate at Antares and IceCube for the 
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two DM scenarios. We pointed out that Antares and IceCube might detect neutrino signals 
in annihilating DM scenario and distinguish different DM scenarios. Section IV contains our 
discussions and conclusions. 



II. NEUTRINO SIGNALS FROM ANNIHILATING/DECAYING DM 

PAMELA data show an excess on positron fraction about 10 ~ 1000 times larger than 
usual estimation, so it is natural to expect that the neutrinos from DM are also enhanced 
by the similar number. In the annihilating DM scenario, it is achieved by Sommerfeld 
enhancement and/or clumps in DM density. In this paper, BF is used to denote the Som- 
merfeld enhancement. Astrophysical enhancement is contained in the astrophysical factor, 
as described below. 

The neutrino flux observed on the Earth can be written as 



and 



>"(E.O)^PefieX^-L-§xA»). (2) 



with dimensionless J (0) and J (0) defined as 



j^{e) = I p\i)di (3) 

PqJ^Q J LOS 

jDi 



J^{9) = -—- / p{l)dl. (4) 
Po-tt© Jlos 

Here 'A' and 'D' denote annihilating and decaying DM respectively, p© = 0.34 GeV cm~^ 
is the local DM density and Rq = 8.5 kpc is the distance between the Sun and the GC. 9 
is defined as the angle between the observational and the GC directions, is the mass of 
DM particle, and dN/dE is the energy spectrum of u per annihilation or decay, (erf) is the 
thermally averaged annihilation cross section, which can be written as 

{av) = {av)Q X BF. (5) 

Here {crv)^ is set to be 3 x 10"^^ cm^s"^, which is the typical annihilation cross section for the 
present dark matter abundance under the standard thermal relic scenario. Thus the number 
BF is enhancement factor compared with the typical value. The BF can arise from the 
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Sommerfeld enhancement as described above. The integral path in Eqs. and (jlj) is along 
the line-of-sight (LOS). The neutrino spectrum dN/dE comes from /i, r decay. We count all 
three flavor of neutrinos and anti- neutrinos at the decay point. After vacuum oscillation, we 



adopt assumption that the three flavor neutrinos have the equal flux . Thus we multiply 
a factor of 1/3 in order to get the flux of muon neutrino and anti-muon neutrino, which is 
the main measured component in the neutrino detector. We should mention here that the 
tau neutrinos can also be observed by the Cherenkov neutrino detector because tau neutrino 
has the different signal compared with muon neutrino. The tau neutrino can interact with 
nucleons and produce tau leptons which subsequently decay into muons. Thus they will 
produce a kink in the reconstruction by Cherenkov light, while the muon neutrino has no 



such kink. The authors of Ref. 



16| have discussed the tau neutrino from decaying dark 



matter and found it can improve the detection by a few orders of magnitude. The reason is 
simply that the atmospheric neutrinos has much less tau neutrinos than muon neutrinos at 
high energy. In this paper, we will concentrate on muon neutrinos observation. 
The solid angle average of this J factor is defined as 



jA,D 

'^^'^ An 



with AQ = 27r(l — cos^). The averaged neutrino flux in a certain cone can be written as 

j,A /pN _ ,2 p ^ 1 {av)dN , 

&iE) = PeRe X X J^^. (8) 

In order to calculate the neutrino flux, we must know both the astrophysical Jaq factor 
and the neutrino spectrum from charged leptons. In Fig. [Tl we plot the J^n factor for 
annihilating and decaying DM respectively. We used NFW |2J] profile here. From the figure 
we can see that the Jaq reach its maximum values, around 1000 and 20 for annihilating and 
decaying DM, in the direction of GC as expected. If cusped density like Moore profile is 
used, one will get much larger and slightly larger J^n factors for annihilating and decaying 
DM respectively. However, some authors 25|] suggested that DM halos around galaxies have 
the shallower density than that of predicted by standard cold dark matter theory. 

Besides the GC, subhalos can also be the high energy neutrino sources. In the direction 
of the DM subhalos, J^n factor can be easily as large as 100, or even higher. In order to 
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FIG. 1: Jaq for annihilating and decaying DM as the function of half-angle for the cone centered 
at the direction of the GC. The DM distribution adopted here is NFW profile. 

keep the line of this paper, we put the analysis on the number distribution of large Jaq 
subhalo in the Galaxy in the Appendix. 

In Fig. [2], the positron fraction from pure muon and tau decays, as well as the energy 
spectrum of positron/electron are plotted for annihilating and decaying DM. It is a model- 
independent discussion of these charged lepton channels, thus we assumed that XX ~^ 
for annihilating DM and x ^ fo^^ decaying DM respectively. From the figures, we can 
see that these channels can account for both the PAMELA and ATIC data with the proper 
parameters. In the following calculation of neutrino flux, we use the same set of parameters 
to get the neutrino spectrum dN/dE. 

Based on the Jaq factor and the neutrino spectrum from charged leptons, we can calculate 
the neutrino flux. In Fig. [3], the expected neutrino flux from muon and tau decay for 
annihilating and decaying DM from the GC are plotted. We choose the half angle 9 to be 
10° which is the angular resolution of Super-K for the neutrino energy from 1 GeV to 10 GeV 



for a conservative analysis [18j. For the neutrino with higher energy, the Super-K is difficult 
to measure the energy of the muons but can count the number. The parameters BF and life 
of DM are taken to be the values which can account for the PAMELA and ATIC data. The 
figures show that the two DM scenarios have comparable average neutrino flux with that of 



the atmospheric neutrino 26|] for high energy neutrino, because of the logarithmic decrease 
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FIG. 2: The positron fraction from pure muon and tau decays, as well as the positron/electron 
energy spectrum for annihilating and decaying DM to account for PAMELA and ATIC data. In 
the figures, the labels on the left, e.g. 0.8 TeV, are the energy of the muon or tau leptons, while 
the labels on the right, e.g. 860 or 0.85, are the boost factor or life in unit of lO^^s for annihilating 
DM and decaying DM respectively. 

of the atmospheric neutrino. For the low energy region, the neutrino flux are much smaller 
than that of the atmospheric neutrino. For some heavy DM, the Super-K can even place 



the constraints [ll, 



n 



Il4| on their annihilating cross section or life [13|. For the smaller cone 
with 6 = 2°, high energy neutrino from the annihilating DM have much larger flux than 
that from the background. 

In the subhalo, the Sommerfeld enhancement for the DM annihilation can be larger than 
that in the halo because the relative velocity of the DM particles are lower. For example 



the subhalos can have the velocity dispersion of O{10)km/s 



which is smaller than the 
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FIG. 3: The muon and anti-muon neutrino flux for annihilating and decaying DM at the GC. The 



atmospheric neutrino flux data is from Ref. 



2a]. The X-axis is the energy of neutrinos and the 



half angle 6 is 10° and 2° respectively. Other labels in the figure are similar with Fig. [2l 

usual velocity of DM v ~ 10~^c ~ SOO/cm/s. Thus one can expect an extra Sommerfeld 
enhancement for the subhalo [6|, 19[ . The larger BF in the subhalo will not contradict with the 
current atmospheric data, since this effect will be diluted with the larger cone angle where 
J factor decreases rapidly. In this paper we do not adopt this extra enhancement, however 
we note here that such enhancement can make the discovery of neutrino from subhalo in 
annihilating DM scenario even more promising. The decaying DM have much smaller J 
factor for massive subhalo, since it does not benefit from the cusped profile very much. 
Thus it is possible to detect neutrino from subhalo in annihilating DM scenario, but very 
challenging in decaying DM scenario. 
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III. MUON RATE AT NEUTRINO TELESCOPES 



The neutrino detector Antares have angular resolution of 2° from lOOGeV to ITeV 27 1 
and can see GC for 63% of a year 2]J. The Antares has about eight years of operation time 
from 2007 to 2015. Using the neutrino effective area A'^^^{E) of Antares 29|], we calculate 
the events in the 2° cone in the Table. [T] for annihilating DM. The total number of events is 
Njj = J A^J^ {E)^jgdE. From the table, we can see that Antares may discover the neutrino 
signal from GC in annihilating DM scenario, but difficult in decaying DM scenario because 
the signal gets negligible increase from smaller cone, which can be seen in the lower panel 
of Fig. [3l The future neutrino detector like KMSNeT can further explore the annihilating 
DM scenario. 
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TABLE I: The neutrino event numbers in the energy interval fyQQGeV — ITeV for eight years of 
Antares operation from the 2° cone in the GC direction, a is the significance defined as S/ v^. 



IceCube locates in the south pole and covers the northern sky. It is excellent to look at 
the possible large DM subhalo in the Galaxy, provided that it has good angular resolution 
like 1° 2^, which can greatly suppress the atmospheric neutrino background. The massive 
DM subhalo can have high density in the center, thus it can produce large neutrino flux and 
can be identified by the high resolution neutrino detector like IceCube. Subhalo can easily 

n 

: 1°) ~ 100 or even larger values p^. In the Appendix, 
the number distribution of massive DM subhalos in our Galaxy has been discussed in detail. 
In the following estimation, we do not adopt the extreme subhalo model, instead we take J 
factor equal to 100. 

For the IceCube, we use the same method in Ref. to estimate the muon rate in 

the neutrino telescopes in annihilating DM scenario. The total muon and anti-muon rate is 



reach large J factor, say J^^'^^°i 
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expressed as following, 



dE„ 



daP{E,^,E^ 
dE„ 



-Pp 



da:{E,^,E^ 
dE„ 



-Pn 



R^{E^)Aeff{E^)dEy^ 



(9) 



where d(j)i,^/dE^^ is the muon neutrino flux arrived at the neutrino detector. 
da^{Ei^^, E^)/dE^ and da^{E^^, E^) / dE^ are differential cross sections for the muon pro- 
duction process vp IX and vn IX. The densities of protons and neutrons near the 
detector are taken to be Pp = |A^^cm~^ and pn = ^NAcm~^ respectively for IceCube where 



detector volume is flUed with ice. Na is the Avagadro's number. The muon energy loss 
is given as ^ = —a — PE, where a and /? are empirical parameters. The distance that 
a muon travels in the Earth before its energy drops below threshold energy E'^'^'^, is called 
muon range which is given as 



(10) 



where we take the par ameters a = 2.0 x lO'^TeVcrn^ / g and P = 4.2 x lO'^cm'^ /g to be 



the same as Ref. 



111]. The effective area of neutrino telescope ^e// is the function of 



muon energy which increases as the energy goes up. The symbol (z/ u) means that we 
also count the anti-muons from the anti-muon neutrinos because the Cherenkov neutrino 
telescope detects both muons and anti-muons. 
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FIG. 4: The total muon and anti-muon rates for annihilating DM in the massive subhalo for the 
IceCube. The half angle 6 is taken to be 1°. 



In Fig. m we show the muon rate as a function of muon energy for several parameters 
which can account for PAMELA and ATIC data. From the flgures we can see that the 
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signal from annihilating DM can be larger than the atmospheric background for heavy DM, 
even without the extra enhancement from lower velocity dispersion. We calculate the muon 
number in the energy interval SOOGeV^ ~ ITeV for ten years and give the significance of 
signal a = S/ \fB in the Tab. [TTl Heavy DM can reach the five sigma significance which 
is quite encouraging for subhalo neutrino search in the annihilating DM scenario. On the 
other hand, if no neutrino hot spots are identified, it can put constraints on subhalo model 
and BF for annihilating DM. 
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TABLE II: The total muon and anti-muon numbers in the energy interval 500GeV^ — ITeV for ten 
years operation of IceCube for massive subhalo. a is the significance defined as 5/ y/B. 

Finally we would like to discuss the limits from the Super-K observation. Super-K have 
searched neutrino signals from DM annihilation in the direction of the Sun, the Earth 
and the GC. Currently, no high energy astrophysical neutrinos have ever been detected. 
Super-K collaboration detected the upward going muons from several GeV to lOTeV in 
the direction of the center with half angles from 5 to 30 degrees, and set limit on muon 
flux l3]. For example, Super-K collaboration placed an upper bound of muon flux in the 
half angle 10° around the direction to the GC as ~ 5 x 10 ""cm 



15^^ 2^1 Q^j, numerical 



examples for annihilating DM, the muon flux are 4.2 x 10 ^^cm "^s ^, 6.5 x 10 ""cm 



-15, 



1.4 X 10 ^^cm 5.1 X 10 ^•'cm' 



-15, 



h-\ 1.7 X 10" 



14 



cm 



and 3.4 x lO'^^ 



cm 



for 



O.STeV /i, ITeV /i, 1.5TeV /i, ITeV r, 2TeV r and 3TeV r respectively. We should emphasize 
that in our calculations we have assumed that DM annihilate into pure muons or tans. If 
DM can annihilate to electron/positrons, it is easy to fit PAMELA and ATIC data with 
much smaller BF for muon and tan channels and avoid the violation of Super-K limits. On 
the other hand, the DM density distribution is still not very clear, especially in the central 
region of Galaxy and subhalos. Different DM profiles give very different contributions to 
<^An(6'), especially when 6 is small. For a cusped profile, the JAn(^) for small 6 will be larger 
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and can compensate the decrease of BF, which can keep the neutrino signal large enough 
and avoid the violation of Super-K limits simultaneously. We would like to mention that 
more stringent constraints on the DM profile may come from the detections of gamma ray 



and synchrotron radiation in the GC 



30|. 



IV. CONCLUSIONS AND DISCUSSIONS 



Based on the PAMELA and ATIC observations, there may exist extra sources of charged 
leptons, namely electron, muon and tau. In this paper, we calculated the neutrino flux 
from muon and tau which arise from the annihilating/decaying DM. The final muon rate at 
neutrino telescopes Antares and IceCube are also estimated. Our results show that Antares 
is promising in discovering the neutrino signal from GC in annihilating DM scenario, but 
challenging in decaying DM scenario. Moreover, massive DM subhalo is also the promising 
high energy neutrino source in the annihilating DM scenario. For ten years operation, 
IceCube can reach five sigma significant neutrino signal from subhalos for heavy DM. Note 
that in the subhalo the extra Sommerfeld enhancement from lower velocity dispersion can 
further improve the neutrino signal. 

In this paper we focus on the muon neutrino detection. Actually the tau neutrino is 
deserved further investigation because the atmospheric neutrino contains less tau neutrino 
than muon neutrino. If one only counts the tau numbers in the neutrino detector, it can 
further suppress the atmospheric background [l^ and this can be done in future kilometer 
neutrino detector like KMSNeT. Finally, it is worthy to mention that pulsars are difficult to 
be detected by upcoming one may distinguish it from the DM 

scenarios. 
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Note added: During the completion of this work one similar analysis appeared 



32|. As a 



cross check, we have calculated the muon flux in the 10 degree half-angle cone from the GC 
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for various channels, which agree with their results. 



APPENDIX A: MASSIVE SUBHALO WITH LARGE Jaq 

In order to determine the subhalo DM contributions to neutrino flux, generally speaking, 
we need to know the DM proflle in the subhalo and the number density of subhalos in the 
Galaxy. Currently N-body simulations provide useful information about the DM distribution 
in the Galaxy. In this Appendix, we utilize the models based on N-body simulations to 
investigate the number distribution of the massive DM subhalo with large J An, which is 
defined in Eqs. and (jS]). 



1. DM subhalo profile 

Based on N-body simulations, the DM distribution can usually be parameterized as, 

P('') = (r/r,)7[l + (r/r,)°](/3-7)/"' ^^^^ 
where ps and are the scale density and scale radius parameters respectively. The param- 



eters {a, 7) are (1, 3, 1) and (1.5, 3, 1.5) for NFW (2^ and Moore [33|] profiles respectively. 

The two free parameters ps and for a subhalo can be determined once we know the 
subhalo mass and the concentration parameter which depends on the specific subhalo 
model. The density scale ps can be determined by using the mass relation 

Ps{r)dV = M,. (A2) 



In the following we will concentrate on how to determine r^. 

If My is known, we can calculate the virial radius r„ of the subhalo which is often ap- 
proximated as the radius within which the average density is greater, by a specific factor 
A = 200, than the critical density of the Universe pc = 139 Mq kpc"^ (M© is mass of the 
Sun). Thus Ty can be expressed as 

r„ describes the radius within which the subhalo can hold together by its own gravity, but 
it does not describe how the DM mass are distributed inside the distance r„. On the other 
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hand, the scale radius describes the DM mass distribution in the subhalo that most of 
mass are within this radius. The concentration models assume that the is proportional 
to the Ty and determine ratio of them. Large ratio of r^/r^ represents the mass of subhalo 
are highly concentrated in a small radius in the center, which results in the cusped profile. 

For a concentration model the ratio of r^/rg is closed related to concentration parameter 
Cy, which is defined as c„ = Here r_2 is another radius parameter which is defined 

as ^(r^p) |r=r_2 = 0. r_2 represents the radius within which the mass of the subhalo is 
concentrated. To illustrate this point, we can easily check that r„2 is in fact the inflexion 



point of the mass, i.e. 



d?M{r) 



r=r_2 



with M(r) = p{r')-A'Kr''^dr' . The relation between 
r„2 and can be calculated by the definition of r„2 and Eq. (]A1|) . The for NFW and 
Moore profiles are r"-^'" = r_2 and r*^°°^^ = r_2/0.63 respectively. Thus the can be 
determined by My, i.e. 

rv{My) ^moore r^^M^) 



.nfw 



(A4) 



c„(M,)' ' 0.63 c„(M„)" 

Based on the above description, we can see that once the Cy — My relation is given by the 
concentration model, the DM profile of subhalo is determined. 



We use the same method as Ref. 
ENSOl 



15 [which adopts two concentration models, which are 



34| and BOl |35|. In the Ref. [36|], the Cy is fitted in a polynomial form as 



In(c^) = y^Cj 



i=0 
-4 



In- 



M 



(A5) 



where Ci = {3.14,-0.018,-4.06 x 10-^0,0} and {4.34,-0.0384,-3.91 x 10-^ -2.2 x 
10-^ -5.5 X 10"^} for ENSOl and BOl model respectively. 

In Fig. 0, the scale radius T g clS db function of subhalo mass My is plotted to show how 
large the subhalos are. The subhalo radius is smaller in BOl model than that in ENSOl 
model because the c„ is larger in BOl model. This plot assumes that subhalos have NFW 
profile. For the Moore profile r™°°''^ is equal to r"-^"'/0.63. 



2. Subhalo number distribution 

The number of subhalos has been revealed by many N-body simulations [37[ . The number 
density of subhalos have a power-law relation with subhalo mass and an uniform distribution 
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FIG. 5: The scale radius as a function of subhalo mass M, 



in the sohd angle as, 



dN 



dM„ ■ Anr'^dr 



No 



M, 



host 



2 ' 



(A6) 



where r is the distance from the GC to subhalo, Mhost is the mass of host halo (the virial 
mass of the Galaxy of Mhost = 10^"^ Mq). th is the core radius which is usually a fraction 
of the virial radius of the host halo (The relation ^ 0.14r^ is adopted here 
slope a varies from 1.7 to 2.1, and an intermediate a is adopted, i.e. a = 1.9 
the normalization factor which is determined by setting the number of subhalos with mass 
larger than lO^M© to be 100 361] . Thus if M„ and the distance r of the subhalo are known, 
we can calculate the number density of the subhalos. 




3. Number of subhalos with JAn{i°) > 100 

In this subsection, we will discuss the number of massive subhalo with JAn{i°) > 100 in 
our Galaxy. From the above discussions, once the concentration model (equivalently saying, 
the function c„(M„) ) is given, the J^n is determined by and the subhalo position in the 
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Galaxy. Generally speaking, the heavier subhalo will have larger Jah- At the same time, 
the distance between the Earth and the subhalo is closer, the Jaq will be larger. 

A monte carlo method is adopted to calculate the number distribution of DM subhalos 
with mass larger than IO^Mq, according to the distribution function Eq. flA6l) . Besides 
My and subhalo distance r from GC, we still need the azimuthal angles 6 and 0, which are 
uniformly distributed in the solid angle to specify the position of the subhalo. Once we 
know the above four parameters, the J^n can be calculated for any specific concentration 
model. The concentration models of ENSOl and BOl are adopted in the calculations, while 
the DM profiles of subhalos used here are the NEW and Moore ones. The total number 
of subhalos N{> IO^Mq) are calculated to be about 6400 in our Galaxy for NEW-ENSOl, 
Moore-ENSOl, NEW-BOX and Moore-BOl subhalo distributions. 

In Eig. El we give the cumulative number of subhalos with JAn{i°) lager than a specified 
value. Erom the figure we can see the number of subhalos with JAn(i°) > 100 depends on the 
concentration model and subhalo DM profile. Moore-BOl model provides the largest number 
of subhalos with high luminosity. There are several subhalos which are Jac(i°) > 100 in our 
galaxy in this model. Eor the Moore-ENSOl and NEW-BOl model, the subhalo number with 
Jaq{i°) > 100 in 6400 subhalos is about 0.6 and 0.1 respectively and even smaller for the 
NEW-ENSOl model. 

In order to show clearly the mass and location of subhalo with JAn{i°) > 100, we depict 
in Eig. [7] the number distribution in the r — My map. Just for illustration purpose, the 
quantity inside the blue contour line shows the total number of these massive subhalos in 
the whole — r map. Outside the blue contour line, the number of such subhalos is zero. 
Moreover the quantity in the red contour line shows the subhalo number in the shadowed 
region. The figure shows that the subhalos tend to have moderate large mass. It is quite 
easy to understand. On one hand, subhalos with larger mass will induce higher J An- On 
the other hand, the number of higher mass subhalos is less. Thus the competition between 
these two effects determines that subhalos with moderate mass tend to provide the largest 
flux. The figure also indicates that subhalos tend to locate around 8.5 kpc from the GC, 
which is the distance between the GC and the Earth. 
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FIG. 6: The cumulative number of sub-halos with astrophysical factor > JAn(i°) a function of 
<^AQ{i°)- The dashed horizon Une corresponds to the case that only one subhalo will be observed. 

4. Summary 

In this Appendix, we showed that the subhalos with JAn{i°) > 100 are very likely to have 
intermediate mass from IO^Mq to IO^^Mq, and locate within several kpc from the Earth. 
The number of these subhalos in our Galaxy is about 0(1) which depends on concentration 
model and subhalo DM profile. It should be noted that the ENSOl and BOl models are for 
distinct halos in the Universe. In our Galaxy, the subhalos are in the host halo (Milky Way 
halo) which is denser than the Universe background. We can expect that the subhalos in 
the Milky Way halo should be more concentrated than the distinct halos. The simulation 
in Ref. 35|] indeed shows the subhalos in a host halo have larger concentration parameter 
Ct, than the distinct ones with the same mass. Thus we can expect larger number of high 
luminosity subhalos in our Galaxy comparing with the numbers above. 

From Fig. [5] and Fig. [71 we can see that the whole subhalo may be difficult to be covered 
by the 1 degree half-cone, which means the subhalo can not be treated as 'point source'. 
Nevertheless, we are interested in the bright core of subhalo which can still be processed as 
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FIG. 7: The number distribution of subhalos with JAr2(i°) > 100- The r is the distance from the 
GC, the My is the subhalo mass. 

point source. The IceCube may discover these hot spots through survey of the sky map by 
1 degree half- cone. 



[1] O. Adriani et al, larXiv:0810.4995] [astro-oh]: [arXiv:0810.4994] [astro-ph]. 
[2] J. Chang et al, Nature 456, 362 (2008). 

[3] D. Hooper, P. Blasi and P. D. Serpico, larXiv:0810.T527l [astro-ph]; H. Yuksel, M. D. Kistler 



19 



and T. Stanev. [a?Xiv:0810.2784l [astro-ph]. 

[4] L. Bergstrom, T. Bringmann and J. Edsjo. [a"rXiv:0808.3725l [astro-ph]: V. Barger, W. Y. Ke- 
ung, D. Marfatia and G. Shaughnessv. larXiv:0809. 01621 [hep-ph]: I. Cholis, L. Goodenough, 
D. Hooper, M. Simet and N. Weiner, larXiv: 0809 1683] [hep-ph]. M. Pospelov and A. Ritz, 
larXiv:0810.1502l [hep-ph]; M. Fairbairn and J. Zupan, ' arXiv:0810.4T47l [hep-ph]. I. Chohs, 
D. P. Finkbeiner, L. Goodenough and N. Weiner, larXiv:0810.5344l [astro-ph]; Y. Nomura 
and J. Thaler, larXiv:0810.5397] [hep-ph]. D. Feldman, Z. Liu and P. Nath, .larXiv:0810.5762l 
[hep-ph]; P. J. Fox and E. Poppitz. [a?Xiv:0811.0399] [hep-ph]. 

[5] M. Cirein and A. Strumia. [arXiv:0808.3867] [astro-ph] : M. Cirelh, M. Kadastik, M. Raidal and 
A. Strumia, arXiv: 0809. 2409) [hep-ph]. 

[6] N. Arkani-Hamed, D. P. Finkbeiner, T. Slatyer and N. Weiner. [arXiv:0810.0713] [hep-ph]. 

[7] P. f. Yin, Q. Yuan, J. Liu, J. Zhang, X. j. Bi and S. h. Zhu, larXiv:0811.0T76] [hep-ph]; 
K. Ishiwata, S. Matsumoto and T. Moroi. [arX"iv:0811.0250l [hep-ph]: C. R. Chen, F. Takahashi 
and T. T. Yanagida, larXiv: 08 1 1 .04771 [hep-ph]; K. Hamaguchi, E. Nakamura, S. Shirai and 
T. T. Yanagida, llarXiv:0811.0737l [hep-ph]; A. Ibarra and D. Tran, larXiv:081 1.15551 [hep-ph]; 
C. R. Chen, F. Takahashi and T. T. Yanagida, larXiv:081 1.33571 [astro-ph]. 

[8] H. Baer, K. m. Cheung and J. F. Gunion, Phys. Rev. D 59, 075002 (1999) 
|arXiv:hep-ph/980636l] ; J. Hisano, S. Matsumoto, M. M. Nojiri and O. Saito, Phys. Rev. D 
71, 063528 (2005) [arXiv:hep-ph/0412403] ; J. Hisano, S. Matsumoto, M. Nagai, O. Saito and 
M. Senami, Phys. Lett. B 646, 34 (2007) [arXiv:hep-ph/0610249 1 . M. Cirelh, A. Strumia and 
M. Tamburini, Nucl. Phys. B 787, 152 (2007) |arXiv:0706.407T1 [hep-ph]]: J. March-Russeh, 
S. M. West, D. Cumberbatch and D. Hooper, JHEP 0807, 058 (2008) |arXiv:0801.3440] [hep- 
ph]]; J. D. March-Russell and S. M. West. larXiv:0812.0559l [astro-ph]. 

[9] M. Lattanzi and J. I. Silk. larXiv:0812.0360l [astro-ph]. 
[10] N. Itoh, H. Hayashi, A. Nishikawa and Y. Kohyama, ApJS 102: 411-424 (1996). 
[11] J. F. Beacom, N. F. Bell and G. D. Mack, Phys. Rev. Lett. 99, 231301 (2007) 
[arXiv:astro-ph/0 608090 1; H. Yuksel, S. Horiuchi, J. F. Beacom and S. Ando, Phys. Rev. 
D 76, 123506 (2007) |arXiv:0707.0T96l [astro-ph]]. 
[12] V. Barger, W. Y. Keung, G. Shaughnessy and A. Tregre, Phys. Rev. D 76, 095008 (2007) 
|arXiv:0708. 13251 [hep-ph]]; V. D. Barger, W. Y. Keung and G. Shaughnessy, Phys. Lett. B 
664, 190 (2008) |arXiv: 0709. 33011 [astro-ph]]. 



20 



S. Palomares-Ruiz, Phys. Lett. B 665, 50 (2008) |arXiv:0712.1937] [astro-ph]]. 
J. Liu, P. f. Yin and S. h. Zhu, Phys. Rev. D 77, 115014 (2008) jarXiv:0803.2T64] [hep-ph]]. 
P. f. Yin, J. Liu, Q. Yuan, X. j. Bi and S. h. Zhu, Phys. Rev. D 78, 065027 (2008) 
|arXiv:0806.3689l [astro-ph]]. 

L. Covi, M. Grefe, A. Ibarra and D. Tran. larXiv:0809.5030l [hep-ph]. 

S. Desai et al. [Super-Kamiokande Collaboration], Phys. Rev. D 70, 083523 (2004) [Erratum- 
ibid. D 70, 109901 (2004)] |arXiv:hep-ex/0404025| . 

Y. Ashie et al. [Super-Kamiokande Collaboration], Phys. Rev. D 71, 112005 (2005) 



|arXiv:hep-ex/0501064|. 



M. Ackermann et al. [The AMANDA Collaboration], Phys. Rev. D 71, 077102 (2005) 



[arXiv : astro-ph /041 2347] ; M. Ackermann et al. [IceCube Collaboration], Astrophys. J. 675, 
1014 (2008) |arXiv:0711.3022] [astro-ph]]. 

M. D. Kistler and J. F. Beacom, Phys. Rev. D 74, 063007 (2006) |arXiv:astro-ph/0607082"] ; 
J. F. Beacom and M. D. Kistler, Phys. Rev. D 75, 083001 (2007) |arXiv:astro-ph/070175l| . 



E. Aslanides et al. [ANTARES Collaboration], arXiv:astro-ph/9907432 



J. Ahrens et al. [IceCube Collaboration], Astropart. Phys. 20, 507 (2004) 
|arXiv :astro-ph/ 0305196] . 

J. Carr et al. [KM3NeT Collaboration], larXiv:0711.2145i [astro-ph]. 

J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 490, 493 (1997) 
|arXiv:astro-ph/9611107] . 



P. Salucci, F. Walter and A. Borriello, arXiv:astro-ph/0206304 



M. Honda, T. Kajita, K. Kasahara, S. Midorikawa and T. Sanuki, Phys. Rev. D 75, 043006 
(2007) [ aFxiv : astro- ph/0611418| . 

T. Stolarczyk [ANTARES Collaboration], Nucl. Phys. Proc. Suppl. 165, 188 (2007). 

L. E. Strigari, S. M. Koushiappas, J. S. Bullock and M. Kaphnghat, Phys. Rev. D 75, 083526 



(2007) |arXiv:astro-ph/0611925| . 

R. Bruijn, PhD thesis. The Antares Neutrino Telescope: Performance studies and analysis of 
first data. 

N. F. Bell and T. D. Jacques, larXiv:0811.082T] [astro-ph]; I. Cholis, G. Dobler, 
D. P. Finkbeiner, L. Goodenough and N. Weiner, I arXiv:081 1.36411 [astro-ph]; G. Bertone, 
M. Cirehi, A. Strumia and M. Taoso. [arXiv:0811.3744] [astro-ph]; E. Nardi, F. Sannino and 



21 



A. Strumia. [a?Xiv:0811.4153] [heo-ph]: J. Zhang, X. J. Bi, J. Liu, S. M. Liu, P. f. Yin, Q. Yuan 
and S. H. Zhu. larXiv:0812.0522l [astro-ph], 
[31] A. Bhadra and R. K. Dey, larXiv:0812.1845] [astro-ph]. 

[32] J. Hisano, M. Kawasaki, K. Kohri and K. Nakavama. larXiv:0812.0219l [hep-ph]. 

[33] B. Moore, T. R. Quinn, F. Governato, J. Stadel and G. Lake, Mon. Not. Roy. Astron. Soc. 
310, 1147 (1999) |arXiv:astro-ph/9903164| . 

[34] V. R. Eke, J. F. Navarro and M. Steinmetz, Astrophys. J. 554, 114 (2001). 

[35] J. S. Bullock et al, Mon. Not. Roy. Astron. Soc. 321, 559 (2001). 

[36] J. Lavalle, Q. Yuan, D. Maurin and X. J. Bi, Astron. Astrophys. 479, 427 (2008). 

[37] G. Tormen, A. Diaferio and D. Syer, Mon. Not. Roy. Astron. Soc. 299, 728 (1998); 
A. A. Klypin, S. Gottlober and A. V. Kravtsov, Astrophys. J. 516, 530 (1999). 
|arXiv:astro-ph/970819l"] ; B. Moore, S. Ghigna, F. Governato, G. Lake, T. Quinn, J. Stadel 
and P. Tozzi, Astrophys. J. 524 (1999) L19; S. Ghigna, B. Moore, F. Governato, G. Lake, 
T. Quinn and J. Stadel, Astrophys. J. 544, 616 (2000). |arXiv:astro-ph/9910166| ; V. Springel, 
S. D. M. White, G. Tormen and G. Kauffmann, Mon. Not. Roy. Astron. Soc. 328, 726 
(2001). [ |arXiv:astro-ph/0012055 |; A. R. Zentner and J. S. Bullock, Astrophys. J. 598, 49 

(2003) . [ arXiv:astro-ph/0304292l ; G. De Lucia et al, Mon. Not. Roy. Astron. Soc. 348, 333 

(2004) . [arXiv:astro-ph/0306205| ; J. Diemand, B. Moore and J. Stadel, Nature 433, 389 (2005) 
|arXiv:astro-ph/0501589] . 

[38] J. Diemand, B. Moore and J. Stadel, Mon. Not. Roy. Astron. Soc. 352, 535 (2004) 

|arXiv:astro- ph/0402160] . 
[39] P. Madau, J. Diemand and M. Kuhlen. [a?Xiv:0802. 22651 [astro-ph] . 



